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This  paper  deals  with  the  control  system  development  for  a  hybrid  energy  storage  system,  consisting  of  a 
battery  and  a  supercapacitor,  for  a  through-the-road-parallel  hybrid  electric  vehicle.  One  of  the  main 
advantages  deriving  from  the  coupling  of  a  battery  and  a  supercapacitor  is  the  possibility  of  reducing 
battery  ageing,  in  addition  to  energy  efficiency  improvements  when  the  system  operates  in  critical 
climate  conditions.  At  the  moment,  no  specific  controller  has  been  proposed  with  the  aim  of  directly 
reducing  battery  wear.  This  paper  presents  a  novel  model  predictive  controller  and  a  dynamic  pro¬ 
gramming  algorithm  including  a  simplified  battery  ageing  model  in  their  formulations.  The  simulation 
results  of  the  model  predictive  controller  and  dynamic  programming  algorithm  are  compared  with  the 
results  deriving  from  a  rule-based  strategy.  The  rule-based  controller  achieves  a  67%  reduction  of  the  root 
mean  square  values  of  battery  current  along  a  selection  of  driving  cycles  in  comparison  with  the  same 
vehicle  equipped  with  battery  only.  In  the  same  conditions  the  battery  peak  current  is  reduced  by  38%. 
The  model  predictive  controller  and  the  dynamic  programming  algorithm  further  reduce  the  root  mean 
square  value  by  6%  and  10%  respectively,  whilst  the  peak  values  are  additionally  decreased  by  17% 
and  45%. 

©  2014  The  Authors.  Published  by  Elsevier  B.V.Onen  access  under  cc  ry-nc-nd  license 


1.  Introduction 

Supercapacitors  are  evaluated  in  automotive  industry  and 
academia  for  hybrid  energy  storage  systems  (HESSs),  in  coopera¬ 
tion  with  batteries  and  DC/DC  converters.  The  typical  properties  [1  ] 
of  Lithium-Ion  batteries  and  supercapacitors  [2]  are  reported  in 
Table  1. 

According  to  the  literature,  the  first  potential  benefit  of  HESSs  is 
represented  by  the  power  loss  reduction  in  the  energy  storage.  In 
fact,  the  energy  efficiency  of  supercapacitors  is  higher  than  for 
batteries  [3,4],  especially  at  significant  currents.  Moreover,  super¬ 
capacitors  allow  regeneration  even  when  the  vehicle  is  working  in 
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critical  ambient  conditions  (i.e.  at  low  temperatures  [5]),  in  which 
the  battery  cannot  operate  in  regenerative  mode.  In  practical  terms, 
from  the  analysis  of  the  average  working  conditions  of  a  vehicle,  the 
supercapacitor  contribution  in  terms  of  enhanced  regeneration  is 
expected  to  be  quite  limited.  Moreover,  in  case  of  HESSs,  DC/DC 
converters  [6]  are  required  to  decouple  the  battery  and  super¬ 
capacitor  voltage  levels.  These  devices  have  an  efficiency  range  that 
at  least  partially  compensates  the  potential  benefit  of  the  super¬ 
capacitor  when  considering  the  whole  system.  Within  the  HESS,  the 
DC/DC  converter  is  usually  the  interface  between  the  supercapacitor 
and  the  DC-link  bus  voltage,  and  manages  the  whole  amount  of 
power  through  the  supercapacitor  [7],  Schupbach  et  al.  [8]  compare 
different  typologies  of  DC/DC  bi-directional  converters,  in  particular 
half-bridge  converters,  single-ended  primary-inductor  converters 
and  Luo  converters,  reaching  the  conclusion  that  the  half-bridge 
layout  represents  the  most  efficient  option.  The  same  conclusion 
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Table  1 

Energy  storage  systems:  typical  properties. 

Battery  Supercapacitor  Unit 

Power  density  ~1  Up  to  20  [kWkg-1] 

Energy  density  150-200  2-10  [Whkg-1] 

Energy  efficiency  75-90  90-98  [%] 


is  reached  by  Ref.  [9],  A  recent  paper  [10]  suggests  the  installation  of 
the  DC/DC  converter  between  the  battery  and  the  supercapacitor, 
and  defines  four  operational  modes  of  the  resulting  system, 
controlled  through  a  rule-based  algorithm.  In  conditions  of  signifi¬ 
cant  power  demand  and  high  vehicle  speed,  the  battery  is  providing 
energy  directly  to  the  motor  inverter.  This  novel  solution  allows  a 
smaller  size  of  the  DC/DC  converter  in  comparison  with  the  con¬ 
ventional  HESS  layout  currently  adopted  in  automotive  industry. 

Ortuzar  et  al.  [11]  report  an  experimental  fuel  consumption 
reduction  of  about  9%  deriving  from  the  adoption  of  a  HESS  within  a 
parallel  hybrid  electric  vehicle,  however  the  specific  prototype  is 
characterized  by  a  lead-acid  battery,  notoriously  less  efficient  than 
modern  Lithium-Ion  batteries.  A  similar  study  [12]  estimates  en¬ 
ergy  consumption  reductions  of  15—20%  through  the  HESS,  but  also 
in  this  case  a  lead-acid  battery  is  considered,  while  in  Ref.  [13]  the 
efficiency  improvement  is  between  4  and  13%  depending  on  the 
driving  schedule.  A  very  recent  paper  [14],  based  on  vehicle  simu¬ 
lations,  shows  significant  energy  consumption  benefits  (from  27% 
to  42%  depending  on  the  driving  schedule)  deriving  from  the 
adoption  of  a  HESS  on  a  series  hybrid  electric  bus.  However,  the  DC/ 
DC  converter  efficiency  does  not  seems  to  be  considered  in  the 
simulation  model  of  the  power  system.  Moreover,  the  implemented 
algorithm  guarantees  the  charge  sustainability  of  the  super¬ 
capacitor  during  the  driving  schedule,  but  no  reference  is  made  to 
the  charge  sustainability  of  the  battery,  nor  to  any  compensation  of 
the  fuel  consumption  for  considering  the  difference  in  battery  state 
of  charge  between  the  initial  and  final  conditions.  Therefore,  a 
rigorous  comparison  of  the  actual  power  losses  for  vehicle  appli¬ 
cations  with  HESS  and  with  battery  only  is  still  required. 

The  second  (but  potentially  more  important  than  the  first  one) 
benefit  of  supercapacitors  is  represented  by  the  reduction  of  bat¬ 
tery  dynamic  loads  (i.e.  currents)  [15],  which  have  direct  impact  on 
battery  life  expectancy.  According  to  the  preliminary  analysis  of 
[16],  based  on  experimental  lead-acid  battery  data  from  Refs.  [17] 
and  [18],  HESSs  bring  a  potential  increase  of  battery  life  between 
137%  (for  traction  batteries)  and  253%  (for  starter  batteries).  How¬ 
ever,  a  direct  estimation  of  modern  Lithium-Ion  battery  life 
extension  due  to  the  adoption  of  a  HESS  is  still  missing  in  the 
literature,  together  with  an  on-line  control  formulation  of  the  po¬ 
wer  split  between  battery  and  supercapacitor,  directly  considering 
battery  wear.  Several  simulation  studies  present  various  forms  of 
empirical  battery  models  [19,20],  which  are  suitable  for  a  compu¬ 
tationally  efficient  estimation  of  battery  wear  for  control  purposes. 
In  general,  battery  wear  depends  on  the  current  profile,  the  depth 
of  discharge  (i.e.  the  same  current  profile  provokes  more  wear  at 
low  state  of  charge)  and  thermal  conditions.  Masih-Tehrani  et  al. 
[14]  adopt  a  specific  wear  model  for  LiFeP04  batteries  in  order  to 
develop  a  dynamic  programming  algorithm  based  on  the  minimi¬ 
zation  of  the  management  cost  of  a  hybrid  electric  vehicle, 
including  fuel  consumption  and  periodic  battery  substitutions. 
However,  the  calendar  life  of  the  battery  is  not  considered  and  the 
amount  of  increase  of  battery  life  is  not  explicitly  presented. 

A  variety  of  controllers  for  HESSs  are  discussed  in  the  literature. 
For  example,  in  the  area  of  heuristic  controllers: 

i)  the  ‘all  or  nothing’  strategy,  which  uses  the  supercapacitor 
(and  not  the  battery)  when  its  state  of  charge  is  above  an 


assigned  critical  level,  and  only  uses  the  battery  when  the 
supercapacitor  is  discharged  [21]; 

ii)  a  strategy  similar  to  i),  however  providing  a  smooth  transi¬ 
tion  between  the  areas  covered  by  one  component  only.  The 
weighting  factor  of  the  power  contribution  of  the  battery  and 
supercapacitor  is  based  on  a  look-up  table  as  a  function  of  the 
supercapacitor  state  of  charge  [21]; 

iii)  a  strategy  based  on  a  bi-dimensional  look-up  table, 
computing  the  power  fraction  to  the  battery  as  a  function  of 
the  states  of  charge  of  the  battery  and  supercapacitor  [21]; 

iv)  the  filtration  strategy,  based  on  the  low-pass  filtering  of  the 
traction  current  for  the  computation  of  the  battery  contri¬ 
bution,  whilst  the  supercapacitor  contribution  derives  from 
the  high-pass  filtering  of  the  traction  current  [21,22]; 

v)  a  vehicle  speed-based  strategy,  according  to  which  the 
supercapacitor  is  charged  up  to  full  capacity  during  periods 
of  low  vehicle  speed  such  that  it  is  ready  to  provide  power 
during  acceleration.  Conversely,  during  periods  of  high 
vehicle  speed,  the  ultracapacitor  should  be  at  the  lowest 
desired  state  of  charge  in  order  to  fully  accept  regenerative 
charge  during  a  braking  event.  An  empirical  formula  estab¬ 
lishing  the  relationship  between  the  supercapacitor  refer¬ 
ence  voltage  and  vehicle  speed  is  provided  in  Ref.  [11].  Ref. 
[23]  adds  a  rule  to  the  speed-based  strategy,  in  order  to 
explicitly  limit  the  peaks  of  battery  current.  Ref.  [24]  uses  the 
same  empirical  formula  for  calculating  the  reference  voltage 
of  the  supercapacitor,  and  then  a  proportional  controller 
based  on  the  difference  between  the  reference  and  actual 
supercapacitor  voltage  is  adopted  for  the  calculation  of  the 
supercapacitor  power; 

vi)  a  rule-based  strategy  based  on  the  combination  of  i)— v),  with 
the  rules  taking  into  account  the  state  of  charge  of  the  two 
components  and  the  HESS  current  demand,  with  an  implicit 
limitation  of  the  battery  current  profile  (despite  no  filter  is 
included).  The  authors  of  this  algorithm  [16]  insist  (without 
quantitative  figures)  on  the  potential  life  extension  of  the 
battery  as  the  main  consequence  of  the  adoption  of  the  HESS 
and  the  proposed  controller. 

In  Ref.  [21  ]  strategies  i) — iii)  are  compared  against  each  other  in 
terms  of  battery  current  profile,  and  iii)  provides  the  smoothest 
profile,  however  the  battery  model  adopted  for  the  assessment  of 
the  controllers  does  not  contain  an  explicit  calculation  of  battery 
wear. 

Borhan  et  al.  [25]  describe  a  model  predictive  controller  (MPC) 
for  the  power  split  between  battery,  supercapacitor  and  internal 
combustion  engine  in  a  hybrid  electric  vehicle.  However,  the  con¬ 
trol  structure  integrates  the  drivetrain  power  split  (i.e.  the  one 
between  engine  and  electric  motor)  and  the  HESS  power  split  into  a 
single  controller,  which  is  not  industrially  viable.  In  fact,  the  same 
hybrid  electric  vehicle  should  undergo  a  redesign  of  its  whole 
powertrain  controller,  when  passing  from  a  battery  only  energy 
storage  to  a  HESS,  without  any  modification  of  the  hybrid  power- 
train.  The  same  limitation  applies  to  the  Pontryagin’s  minimum 
principle-based  controller  presented  in  Ref.  [12],  which  is  run  off¬ 
line  for  computational  reasons.  Moreover,  the  MPC  in  Ref.  [25]  is 
designed  for  tracking  a  reference  state  of  charge  of  the  super¬ 
capacitor,  which  is  not  clearly  specified  in  the  paper.  Finally,  Ref. 
[11]  proposes  a  neural  network-based  strategy,  with  the  specific 
purpose  of  reducing  the  power  losses  in  the  HESS. 

The  conclusion  is  that  there  is  no  wide  consensus  regarding  the 
assessment  of  the  actual  main  benefit  deriving  from  the  adoption  of 
the  HESS  (efficiency  or  battery  wear  reduction),  nor  any  on-line 
HESS  management  algorithm  designed  to  explicitly  take  into  ac¬ 
count  battery  wear. 
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The  objectives  of  this  paper  are: 

•  To  formulate  an  innovative  model  predictive  controller  for 
hybrid  electric  or  fully  electric  vehicle  HESS  management,  spe¬ 
cifically  designed  for  the  reduction  of  battery  wear; 

•  To  design  a  novel  dynamic  programming  (DP)  formulation  for 
estimating  the  ideally  achievable  increase  of  battery  life  dura¬ 
tion  through  the  HESS.  Simplified  forms  of  DP  can  become  viable 
in  future  by  using  the  information  from  the  cloud  [26]  and  the 
vehicle  navigation  system  for  programming  the  schedule  of  the 
whole  trip; 

•  To  assess  the  performance  benefits  (i.e.  battery  life  extension 
and  power  losses)  of  MPC  and  DP  against  a  computationally 
efficient  rule-based  algorithm  from  the  literature  (formulation 
vi)  discussed  in  this  section),  and  outline  the  trade-off  between 
model-based  controllers  and  heuristic  formulations  for  HESS 
management; 

•  To  compare  the  energy  storage  power  losses  in  case  of  battery 
only  and  of  HESS,  for  a  case  study  through-the-road-parallel 
(TTRP)  hybrid  electric  vehicle  (HEV)  along  a  selection  of 
driving  cycles. 

The  paper  presents  the  simulation  model  including  the  high- 
level  controller  for  the  power  split  between  battery  and  super¬ 
capacitor,  three  formulations  of  HESS  management  controllers  for 
the  increase  of  battery  life  expectancy,  and  a  detailed  analysis  of  the 
simulation  results. 


Table  2 

Main  vehicle  parameters. 


Symbol  Description 


Quantity  Unit 


M  Overall  vehicle  mass  (unladen)  2036 

1  Wheelbase  2.99 

hcG  Height  of  the  centre  of  mass  0.47 

Ca  Aerodynamic  drag  coefficient  0.35 

S  Frontal  area  2.24 

«EMmax  Maximum  electric  motor  speed  12,000 

TEMmax  Maximum  electric  motor  torque  120 


McEmax  Maximum  internal  combustion  engine  speed  8000 

hcEmax  Maximum  internal  combustion  engine  torque  683 


[kg] 

[m] 

[m] 

1-1 

[m2] 

[rpm] 

[Nm] 

[rpm] 

[Nm] 


TWice  is  given  by: 

TWkeM  =  TlCEtGB,CEtDiffICEJ?GB,CETlDiffIcE  ~  Je qICE^W  (2) 

where  Tice  is  the  engine  torque,  tGBice is  the  ICE  axle  gearbox  ratio, 
TDiff,CE 's  the  ICE  differential  gear  ratio,  i/gb1ce  >s  the  ICE  gearbox  ef¬ 
ficiency,  iJoifficE  is  the  ICE  final  reduction  efficiency,  JeqICF  is  the 
equivalent  inertia  of  the  ICE  drivetrain  at  the  wheel  and  1%  is 
wheel  acceleration.  A  first  approximation  of  Je qiCE  can  be  calculated 


JeqICE  =  Jw  +JDiff,cE  +isS,cETDiffICE??Diff,CE 
+ JPSICE  TGBICE  TDiffICE  7?GBICE  ^Diffe 


(3) 


2.  Simulation  model  and  high-level  controller 

2.1.  Vehicle  model 


where  Jw  is  the  wheel  inertia,  JDiff1CE  is  the  ICE  differential  assembly 
inertia,  Jss]CE  is  the  ICE  secondary  shaft  equivalent  inertia  and  JPS|CE  is 
the  ICE  primary  shaft  equivalent  inertia.  Similarly,  TwEM  is  given  by: 


The  case  study  vehicle  is  a  high  performance  passenger  car 
consisting  of  a  front  electric  axle  equipped  with  a  2-speed  trans¬ 
mission  system  [27]  and  a  rear  engine-driven  axle  equipped  with  a 
7-speed  dual  clutch  gearbox  [28],  The  vehicle  layout  is  shown  in 
Fig.  1,  while  the  main  vehicle  parameters  are  reported  in  Table  2. 
The  simulator  (Fig.  2)  is  based  on  a  hybrid  backward/forward-facing 
approach  [29]  that  allows  the  calculation  of  the  amount  of  energy/ 
power  required  from  the  HESS  along  the  driving  schedule. 

The  overall  drivetrain  torque  at  the  wheel,  Tw,  can  be  calculated 
as  the  sum  of  the  internal  combustion  engine  (ICE)  drivetrain  tor¬ 
que,  TW|CE,  and  the  electric  motor  (EM)  drivetrain  torque,  TWem: 


TWm(t)  =  7'EM'rGBEM TDiff em ''IgBem ??DiffEM  ~JeqEM#W  (4) 

where  Tem  is  the  electric  motor  torque,  tGBem  is  the  EM  axle  gearbox 
ratio,  TDiffEM  is  the  EM  differential  gear  ratio,  7]GBem  is  the  EM  gearbox 
efficiency,  ?jDiffEM  is  the  EM  final  reduction  efficiency  and  JeqEM  is  the 
equivalent  inertia  of  the  EM  drivetrain  at  the  wheel.  A  first 
approximation  of  JeqEM  can  be  calculated  as: 

JeqtM  =  Jw  +.foiffEM  +  JsSEMtDiffEM2j?DiffEM 

+ JPSEM  tGBem  2  TDiffEM  2  ’IgBjm  4DiffEM  (  5  ) 


Tw  (t)  =  Twm{t)  +  TWia(t) 


(1) 


where  JDiffEM  is  the  EM  differential  inertia,  JSsEM  is  the  EM  secondary 
shaft  equivalent  inertia  and  JPSem  is  the  EM  primary  shaft  equivalent 
inertia.  The  drivetrain  efficiencies  are  implemented  in  the  form  of 
look-up  tables  (experimentally  derived  by  the  respective  manu¬ 
facturers)  as  functions  of  input  torque,  speed  and  operating  tem¬ 
perature.  The  wheel  and  vehicle  balance  equations  are  omitted  for 
brevity.  Gearshift  maps  are  adopted  for  the  transmission  gear 
selection. 

2.2.  High-level  controller 

The  high-level  controller  has  to  determine  the  power  split  be¬ 
tween  the  ICE  driven  axle  and  the  EM  driven  axle  of  the  TTRP  HEV. 
The  adopted  algorithm  is  the  equivalent  consumption  minimiza¬ 
tion  strategy  (ECMS)  presented  in  Refs.  [30,31],  Alternative  for¬ 
mulations  are  presented  and  discussed  in  the  literature  [32,33], 

The  core  of  the  strategy  is  the  cost  function  (Eq.  (6)),  defined  as 
the  sum  of  the  fuel  cost  (dimensionally  a  fuel  mass  flow  rate)  and 
the  equivalent  cost  of  the  electric  energy.  At  each  time  step  the 
costs  for  all  the  feasible  splits  are  evaluated  and  the  power  split 
with  the  minimum  cost  is  selected. 


Fig.  1.  TTRP  HEV  layout. 
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Fig.  2.  Schematic  of  the  vehicle  simulator. 


C(t)  =  CEM(t)  +  C,CE(t)  (6) 

where  CEm(0  (Eq.  (7),  in  which  the  two  formulations  refer  to  pos¬ 
itive  and  negative  motor  torques)  is  the  cost  related  to  the  electric 
drivetrain  and  QceM  (Eq.  (8))  is  the  cost  related  to  the  ICE 
drivetrain: 

fSFC„PEM(rEM(t),WEM(t)) 

Hi  2  _  (7) 

SFCRPEM(TEM(t),  wEm(C))??hEsst?em 

CICE(t)  =  BSFCMAP(TICE([),0JICE(t))-P,CE(t)  (8) 

where  5)Hess  and  are  respectively  the  lumped  average  efficiency 
of  the  HESS  and  EM  drivetrains,  SFCr  and  SFCd  are  respectively  the 
equivalent  cost  factors  during  HESS  charge  and  discharge,  PEm  is  the 
EM  power,  Tice,  mice  and  Pice  are  respectively  the  ICE  torque,  speed 
and  power,  and  BSFCMap  is  the  ICE  brake  specific  fuel  consumption. 

The  algorithm  selects  the  feasible  split  with  the  minimum 
cost,  C: 

C*  =  min[C(t)]  (9) 


has  already  been  experimentally  validated  for  several  commercial 
Lithium-Ion  batteries  (e.g.  Sony  US18650,  Panasonic  CGR18650) 
and  because  of  its  low  computational  cost. 

The  model  is  based  on:  i)  an  equilibrium  potential,  EBatt,  defined 
by  a  look-up  table  as  a  function  of  temperature  and  state  of  charge; 
ii)  two  resistors,  Pi.Batt  and  R2,Batt,  for  the  computation  of  the  power 
losses:  and  iii)  a  capacitance,  CBatt,  which  characterizes  the  tran¬ 
sient  response  of  the  battery.  The  state  of  charge  (SOCfiatt)  is 
computed  as: 

SOCgatt  =  Q  f  </Batt)0(TBatt)Wt)dt  +  SOC0,  Batt  (10) 

0 

where  Q  is  the  battery  capacity,  a(/Batt)  takes  into  account  the 
charge/discharge  rate  effect,  /3(TBatt)  takes  into  account  the  thermal 
effect,  fBatt  and  TBatt  are  respectively  battery  current  and  tempera¬ 
ture,  and  SOCo.Batt  is  the  initial  state  of  charge.  The  main  battery  cell 
parameters  (Dow  Kokam  XALT  8  Ah)  considered  in  this  activity  are 
reported  in  Table  3. 

Since  the  battery  behaviour  is  temperature  dependent,  a  ther¬ 
mal  model  is  required.  An  approximated  thermal  energy  balance 
has  been  adopted  [35], 


2.3.  Hybrid  energy  storage  system 

The  HESS  model  includes  the  battery  and  supercapacitor 
models,  which  are  connected  in  parallel  with  the  DC/DC  converter, 
as  shown  in  Fig.  3. 

2.3.1.  Battery  model 

The  battery  model  (Fig.  4)  implemented  in  the  vehicle  simulator 
is  described  in  Ref.  [34],  This  model  has  been  selected  because  it 


-  =lBatt(f)  '^2,  Batt  + 


1 


Rl,  Batt 

[^Batt(t)  -  EBatt(l)  -  R 2,  BattlBatt]2 
-  MPBatt(t)  -  7a] 


(11) 


where  m  is  the  mass  of  the  battery,  cP  is  the  specific  heat  ca¬ 
pacity,  TBatt  is  the  battery  temperature,  VBatt  is  the  potential 


Fig.  3.  Simulated  HESS  layout. 


Fig.  4.  Equivalent  battery  model. 
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Table  3 

Battery  cell  parameters. 


Description  Quantity  Unit 


Nominal  cell  capacity  8  [Ah] 

Nominal  cell  voltage  3.7  [V] 

Peak  cell  discharging  current  24  [A] 

Number  of  cells  in  parallel  6  [-] 

Number  of  cells  in  series  82  J 


across  the  battery  terminals,  hc  is  the  convective  heat  transfer 
coefficient,  A  is  the  equivalent  surface  and  Ta  is  the  ambient 
temperature. 

To  evaluate  the  battery  wear  during  the  driving  cycle,  a  damage 
accumulation  model,  namely  “accumulated  Ah-throughput”  [36], 
has  been  implemented.  It  is  based  on  the  total  amount  of  charge 
that  can  flow  through  the  battery  before  it  reaches  the  end  of  its 
life.  The  model  relies  on  the  definition  of  a  severity  factor  a  (Eq. 
(12)),  which  depends  on  /Batt,  TBat t  and  SOCBatt: 

/•EOL 

,  7(WW  SOCBatt)  Jo  |,Batt|df 

- f  - 7EOL 

Jo  |/n°m|dt 

where  7(/Batt,TBatt,SOCBatt)  is  the  battery  duration  (Ah-throughput) 
corresponding  to  its  actual  operating  conditions  in  terms  of  current, 
temperature  and  state  of  charge,  T  is  the  total  Ah-throughput  when 
the  battery  undergoes  nominal  cycles,  as  defined  by  the  manufac¬ 
turer  (in  the  case  study  the  nominal  conditions  are  1-C  rate,  25  °C 
and  100%  depth  of  discharge).  /n0m  is  the  nominal  current  and  EOL 
indicates  the  end  of  life.  The  severity  factor  represents  the  relative 
ageing  effect  with  respect  to  the  nominal  cycle,  and  is  higher  than  1 
for  conditions  which  are  more  demanding  in  terms  of  battery  wear. 
A  qualitative  example  of  severity  factor  map  is  provided  in  Fig.  5. 
Low  operating  temperatures,  low  values  of  SOCBatt  and  high  C-rates 
provoke  a  significant  increase  of  the  severity  factor.  These  maps  can 
be  obtained  through  specific  experimental  procedures,  which, 
however,  are  beyond  the  common  practices  of  battery  manufac¬ 
turers  [36],  This  uncertainty  in  the  severity  factor  determination  is 
the  reason  why  the  benefits  deriving  from  the  adoption  of  a  HESS 
are  usually  presented  in  terms  of  reduction  of  battery  current,  and 
not  directly  in  terms  of  battery  life.  However,  the  analysis  of  the 
HESS  performance  should  be  actually  based  on  the  extension  of 
battery  life. 


Fig.  6.  Equivalent  supercapacitor  circuit 


To  evaluate  the  actual  depletion  of  the  battery  charge,  the 
following  equation  is  adopted: 

Aheff  =  jj  ^(^Batti  ^Batt-  SOCBatt ) I ^Batt (T) | dT  (13) 

0 

which  represents  the  amount  of  charge  that  would  need  to  go 
through  the  battery  using  the  nominal  cycle  to  have  the  same 
ageing  effect  of  the  actual  conditions.  The  end  of  life  is  reached 
when  Aheff  =  r  [36], 

2.3.2.  Supercapacitor  model 

The  supercapacitor  is  modelled  by  a  series  resistance  RESR.scap, 
the  main  responsible  for  the  losses  in  the  component,  a  capacitance 
Cscap,  and  a  resistor  J?EPR,Scap  to  model  the  self-discharge  current. 
The  capacitor  is  characterized  by  a  variable  capacitance  which  is 
dependent  on  its  cross  voltage. 

This  model  (Fig.  6)  has  been  chosen  because  of  its  flexibility  and 
low  computational  cost,  nevertheless,  it  is  well  known  in  literature 
[37,38].  The  parameters  of  the  circuit  (referring  to  the  NESSCAP 
supercapacitor  5000  F),  included  in  Table  4,  can  be  derived  from  a 
set  of  charge-discharge  tests  [39], 

2.3.3.  DC/DC  converter  model 

In  order  to  take  in  account  the  losses  in  the  DC/DC  converter  an 
efficiency  map  has  been  adopted.  This  map  is  a  function  of  the 
supercapacitor  voltage  and  power  demand  [22], 

3.  HESS  controllers 

3.1.  Rule-based  strategy  [16] 

The  flow-chart  describing  the  rule-based  strategy,  presented  in 
Ref.  [16],  is  reported  in  Fig.  7.  The  algorithm  decides  the  power  split 
depending  on  the  overall  power  request  and  the  state  of  charge  of 


Fig.  5.  Battery  severity  factor  map. 
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Table  4 

Supercapacitor  cell  parameters. 

Description  Quantity  Unit 

Nominal  cell  capacitance  5000  [F] 

Rated  cell  voltage  2.7  [V] 

Peak  cell  current  2547  [A] 

Number  of  cells  in  series  112  [— ] 


the  supercapacitor.  Pmjn  is  the  battery  power  threshold  (in 
discharge)  below  which  the  whole  power  demand  is  met  by  the 
battery,  Pch  is  the  charging  power  sent  from  the  battery  to  the 
supercapacitor,  and  l/req  is  the  required  supercapacitor  voltage, 
which  is  a  decreasing  function  of  vehicle  speed. 

Rule-based  controllers  are  easily  implementable;  however  sig¬ 
nificant  care  is  required  in  the  empirical  tuning  of  their  parameters. 
Some  examples  of  sensitivity  analysis  of  the  impact  of  the  main 
tuning  factors  are  discussed  here.  The  increase  of  predicted  battery 
life  as  a  function  of  Pmin  is  shown  in  Fig.  8.  The  battery  life  peaks  at 
around  6000  W  for  all  the  simulated  driving  cycles,  value  at  which 
the  battery  root  mean  square  (RMS)  current  also  shows  a  minimum 
(Fig.  9).  Finally,  the  maximum  value  of  the  battery  current  for  all  the 
driving  schedules  simulated  as  function  of  Pmjn  is  reported  in 
Fig.  10.  In  summary,  a  good  trade-off  for  the  parameter  Pmin  is 
represented  by  6000  W,  therefore  this  value  has  been  selected  for 
the  final  simulations  and  comparison  with  the  other  controllers. 
Moreover,  a  sensitivity  analysis  of  PCh  has  been  carried  out  (here 
omitted  for  brevity).  Flowever,  its  impact  on  the  predicted  battery 
life,  RMS  and  maximum  current  can  be  neglected. 

In  summary,  the  main  controller  parameters  adopted  in  this 
study  are  reported  in  Table  5.  The  sampling  time  of  the  controller 
has  been  set  at  the  same  value  as  the  sampling  time  of  the  HESS 
model  (0.01  s). 
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Fig.  8.  Battery  life  increase  as  function  of  the  tuning  parameter  Pmin  for  the  different 
driving  cycles. 


makes  the  controller  more  practical  for  a  real  world  vehicle 
implementation  than  the  one  presented  in  Ref.  [25],  including  the 
overall  energy  management  and  the  HESS  management  in  a  single 
control  structure. 

3.2.1.  Battery  model 

The  battery  model  for  MPC  design  is  reported  in  Fig.  12.  It  is 
based  on  the  battery  equilibrium  potential,  Efiatt.  and  an  equivalent 
internal  resistance,  RBatt-  The  component  dynamics  can  be  derived 
from  the  battery  power,  PBatt: 


^Batt  —  ^Batt^Batt  +  ^Batt^Batt  (14) 


3.2.  Model  predictive  control 

The  MPC  algorithm  involves  three  steps  (Fig.  11):  i)  prediction  of 
the  future  outputs  over  the  optimisation  horizon,  using  a  simplified 
model  of  the  system;  ii)  evaluation  of  the  cost  function  for  the  set  of 
future  outputs  of  the  system;  and  iii)  adoption  of  the  control  policy 
with  the  minimum  cost  and  compliant  with  the  system  constraints. 

The  controller  relies  on  a  model  of  the  system  and  a  cost  func¬ 
tion,  which  represent  the  core  of  the  algorithm.  The  model  has  to 
predict  the  system  response  with  good  accuracy,  while  keeping  a 
low  computational  cost.  In  this  study,  reduced  order  models  of  the 
HESS  components  have  been  adopted  as  reported  in  the  following 
paragraphs. 

The  proposed  MPC  has  a  lower  hierarchical  level  than  the  ECMS 
controller  in  the  overall  system  of  Fig.  2.  This  hierarchical  structure 


-£Batt  +  \j ^Batt  +  4RBattPBatt 


_  ^Batt  _ 

Q 


+  y^Lt 


3.2.2.  Supercapacitor  model 

The  supercapacitor  dynamics,  deriving  from  the  simplified 
model  in  Fig.  13,  can  be  expressed  as: 

Pscap  =  ^Scaplscap  +  ^C^Scap  (17) 


Fig.  7.  Rule-based  strategy  flow-chart  [16], 
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Table  5 

Rule-based  controller  parameters. 

Symbol  Quantity  Value  Unit 

Pch  Minimum  battery  power  800  [W] 

Pmln  Power  delivered  by  the  battery  to  the  supercapacitor  6000  [W] 

SOCmax  Maximum  supercapacitor  state  of  charge  1  [-] 

SOCmin  Minimum  supercapacitor  state  of  charge  0.5  [-] 


At  each  sample  time  k  the  system  equations  and  the  charac¬ 
teristics  (i.e.  EBatt(SOCBatt))  of  the  components  are  linearized  around 
the  operating  conditions  at  the  beginning  of  the  control  window. 
The  dynamics  of  the  system  can  be  expressed  as: 


Fig.  9.  Root  mean  square  battery  current  as  function  of  the  tuning  parameter  Pmin  for 
the  different  driving  cycles. 


— SOCscapVc, 
*Scap  =  - 


V(SOCScapVc, 

2^Scap 


max)  +  4RscapPscap 


(18) 


f  SOCscap  1 
l  SOCeatt  J 


^Scap  0  1  f  SOCScap 
o  r^BattJ  1  SOCeatt 


(21) 


while  the  outputs  of  the  system  are  formulated  as: 


dt 

max  +  V (SOCscapUc  n 


f  SOCscap  1 
1  O^Batt  J 


Cscap  0  1  /  SOCscap 
0  Ceatt  J 1  SOCeatt 


0  0  1 
0  Deatt  J 


(22) 


where  Kscap.  ^Batt  and  *Ltt  are  constant  matrices  deriving  from 
the  system  linearization.  By  defining  the  matrices  in  Eq.  (23)  and  in 
Eq.  (24): 


where  fscap  is  the  supercapacitor  current,  Vq  is  the  voltage  across 
the  capacitance,  Rscap  is  the  supercapacitor  resistance,  SOCscap  is  the 
supercapacitor  state  of  charge  and  Vc  max  is  the  maximum  capaci¬ 
tance  voltage. 

3.2.3.  Controller  implementation 

Within  the  implicit  MPC  algorithm  a  state-space  representation 
of  the  system  has  been  implemented  [40],  The  states  x,  the  inputs  u 
and  the  outputs  y  of  the  state-space  formulation  are  reported  in 
Eq.  (20). 

'[-{s«S}'u-{P£}J'-{sSt}  (20) 


B 


E  a‘b 

At„ 

i  =  0 

,  T  = 

£  A'B 

Ah 

D 

E  A'B 

L  i  =  0  J 

(23) 


(24) 


Fig.  10.  Maximum  battery  current  as  function  of  the  tuning  parameter  Pmi„  for  the 
different  driving  cycles. 


where  tw  is  the  length  of  the  control  window  and  tp  is  the  predic¬ 
tion  time,  the  free  evolution  of  the  system  states  over  the  prediction 
time,  x(k),  can  be  expressed  as: 

x(k)  =  k[W{x }  +  rpBatt(k  -  1)  +  K]  +  Xo  (25) 

xo  contains  the  initial  conditions  at  the  beginning  of  the  control 
window,  k  is  the  time  instant  over  the  prediction  time  and  K  con¬ 
tains  constants  deriving  from  the  linearization  of  the  system. 


Fig.  11.  MPC  algorithm  layout 
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attery  model  for  MPC  controller  design. 
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where  j  indicates  a  sample  of  the  prediction,  which  is  discretised  at 
time  intervals  At  (they  can  be  different  from  the  discretisation  time 
adopted  for  k),  ts  represents  the  beginning  of  the  control  window, 
Y(j )  is  the  actual  output  of  the  system,  T(j )  is  the  reference  output  of 
the  system,  and  Q  and  R  are  two  weighting  matrices. 

By  defining  e(k): 


e(k)  =  T(k)  - 


ab(k) 

i 

,<hr(k). 


[Ax(k)+SPBatt(k-l)] 


(27) 


where  <rb  is  the  battery  severity  factor  (Eq.  (12))  evaluated  for  the 
free  response  of  the  system,  the  second  term  on  the  right-hand  side 
of  Eq.  (27)  represents  the  free  response  of  the  system  outputs.  The 
matrices  A  and  S  are  reported  in  Eq.  (28): 


A  = 


0 


0 


°1 
0  ’ 
cj 


(28) 


rr m  is  the  average  value  of  the  severity  fa 


the  prediction 


Jmp c(k)  =  e{k)Tde{k)  -  APBatt(k)TG  +  APBatt(k)T  LAPBatt(k)  (32) 

where  APBatt  is  the  battery  power  difference  between  two 
consecutive  time  samples.  G  is  defined  as  [40]: 

G  =  2  6TQ_e(k)  (33) 

and  L  as: 

L=6TQ_d  +  R  (34) 

Therefore,  the  cost  function  (Eq.  (32))  is  the  sum  of  the  free 
evolution  of  the  system  (e(k)TQe(k))  and  the  term  to  be  minimized 
(APBM(k)TG  +  APBatt(fc)TfAPBatt(k)),  physically  representing  the 
forced  response  of  the  system.  This  results  in  a  quadratic  pro¬ 
gramming  problem  (Eq.  (35))  for  which  efficient  on-line  compu¬ 
tational  techniques  are  available,  such  as  interior  point  [41  ]  and 
active  set  [42]: 


where  the  matrices  C  and  D  (from  the  state-space  formulation  of  ,* 

Eq.  (22))  are  expressed  in  Eq.  (29):  JmpC’  0ptl  ’ 


min  [lAPBatt(k)THAPBatt(k)  +/TAPBatt(k)]  (35) 


(29) 


The  forced  response  E  of  the  system  outputs  can  be  written  as: 


E  =  6APBatt(k) 


(30) 


Aconstr^Wk)  <  Bconstr  (36) 

where  Aconstr  and  BConstr  are  the  two  constraint  matrices  (i.e. 
defining  the  minimum  and  maximum  state  of  charge  for  the 
supercapacitor  and  the  battery),  and  H  and  /  are  two  matrices 
deriving  from  the  quadratic  programming  formulation  [40], 

The  reference  values  for  the  controller  in  Y(j)  (Eq.  (26))  are  the 
supercapacitor  state  of  charge  (SOCscap)  at  the  beginning  of  the 
control  window,  to  provide  a  smooth  power  profile  of  the  super¬ 
capacitor,  and  zero  “corrected"  battery  current,  to  maintain  (r/Batt  as 
small  as  possible.  The  physical  principle  of  this  novel  controller  is  to 
reduce  battery  wear  through  a  reduction  of  the  “corrected”  battery 
current.  The  value  of  the  battery  weighting  factor  in  Q  is  equal  to  1. 
The  supercapacitor  weighting  factor  is  a  function  of  SOCscap.  to  limit 
and/or  prevent  a  quick  discharge  of  the  supercapacitor.  For  the 
simulations  presented  in  this  paper  the  weighting  factor  matrix  R 
has  been  set  at  zero. 

Once  the  optimal  APBatt  has  been  calculated,  the  algorithm 
computes  the  optimal  battery  power: 

fWopt(k)  =  APBatt_0pt(k)+PBatt(k-l) 


(37) 
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Fig.  14.  Battery  life  increase  as  function  of  the  MPC  prediction  time  for  the  different 
driving  cycles. 

The  sample  time  of  the  MPC  has  been  set  at  0.5  s  in  order  to 
achieve  a  low  computational  cost,  comparable  to  the  one  of  the 
rule-based  algorithm.  During  the  controller  implementation  phase, 
a  comparison  of  the  same  MPC  for  the  sample  times  of  0.5  s  and 
0.01  s  has  shown  a  negligible  variation  of  the  results,  with  a  sig¬ 
nificant  increase  of  the  required  simulation  time  in  the  case  of 
0.01  s.  Moreover,  the  battery  life  increase  predicted  for  the  different 
driving  missions  as  a  function  of  the  prediction  time  and  control 
horizon  is  reported  in  the  sensitivity  analyses  of  Figs.  14  and  15.  The 
predicted  battery  life  peaks  for  prediction  times  of  about  10  s  for  all 
the  driving  schedules,  while  it  experiences  an  asymptotic  behav¬ 
iour  as  function  of  the  control  horizon.  Following  this  analysis,  the 
prediction  time  and  the  control  horizons  have  been  set  at  10  s. 

3.3.  Dynamic  programming 

A  dynamic  programming  algorithm  (Fig.  16)  has  been  imple¬ 
mented  to  understand  the  potential  benefits  of  a  global  optimiza¬ 
tion  through  the  whole  driving  cycle.  The  discretization  has  been 
carried  out  according  to  the  supercapacitor  state  of  charge,  SOCscap- 
The  costs  at  each  grid  point  are  calculated  for  each  discrete  time  (k) 
and  consequently  the  path  with  the  minimum  cost  (i.e.  optimal 
control  policy)  is  selected  proceeding  backwards  [43], 

The  cost  function  (Eq.  (38))  is  defined  to  minimize  the  overall 
battery  wear  along  the  driving  cycle: 
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Fig.  16.  Dynamii 

c  programming  schematic  with  a  simplified  grid  (k 
:ates  the  value  of  SOCScap). 

indicates  the  time 

to  evaluate  the  cost  for  each  discretization  point  (in  terms  of  time 
and  SOCscap),  the  supercapacitor  power  is  calculated  according  to: 

Pscap (k;  m.  i) 

—  Pscap_cap(k;  i)  +  Fscap_res(k,  tn, !) 

(39) 

where  Pscap_cap(k,m,i),  defined  in  Eq.  (40),  is  the  power  at  the 
capacitance  required  for  the  variation  of  SOCscap  from  the  value 
corresponding  to  the  generic  discretization  point  i  at  the  sample 
time  k  -  1  to  the  discretization  point  m  at  the  sample  time  k. 
Pscap_res( k,m,i),  Eq.  (41 ),  is  the  consequent  power  dissipation  on  the 
supercapacitor  resistor  due  to  the  variation  of  SOCscap- 

f5Scap_cap(k,  tn,  i)  —  2  ^Scap  |^(SOCscap(k.  m)Vc.  max) 

-  (SOCScap(k-l,l)Vc.  max)2J 

(40) 

^Scap_res(^i  —  ^Scap  [Qcap  (SOCscap (^j  mWc,  max 

~~~  /f.  ^  .vW  m2 

(41) 

-SOCScap(k-l,0W,max)]2 


Subsequently,  the  battery  power  can  be  expressed  as: 


JDP,  optimal  —  tnin  |  if/BaU:  ^Batt-  SOCeattil^battWlj  (38) 

where  <7(/Batt.TBatt.SOCBatt)  is  the  severity  factor  (Eq.  (12)). 

The  models  adopted  for  the  DP  control  action  are  the  same  as 
the  simplified  ones  adopted  for  the  MPC  (Figs.  12  and  13).  In  order 
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Fig.  15.  Battery  life  increase  as  function  of  the  MPC  control  horizon  for  the  different 
driving  cycles. 


Pb att(fc,  rn,  i)  =  PhessW  -  PscaP(k,  rn,  i)y  (42) 

y  =  Vdc/dc  (43) 

where  jjdc/dc  is  the  lumped  efficiency  of  the  DC/DC  converter,  and 
r  =  ±1  depending  on  the  sign  of  the  supercapacitor  power  Pscap-  The 
power  profile  of  the  hybrid  energy  storage  PhessM  is  derived  from  a 
simulation  with  the  vehicle  model  and  the  ECMS  high-level 
controller  detailed  in  Section  2.  Once  the  battery  power  has  been 
determined,  the  current  (Eq.  (44))  and  the  severity  factor  (Eq.  (45)) 
can  be  evaluated  in  order  to  calculate  the  equivalent  cost. 


h att(fc,  ™,  0 


-PBatt  +  \f  Pfiatt  +  4PBattPBatt(ic  Ft,  i) 
2PBatt 


(44) 


o  =  er(katt(k,  m,  i),  TBatt(k,  m,  t),  SOCBatt(k,  m,  i))  (45) 

The  number  of  discrete  points  (range  of  m)  for  SOCscap  has  been 
set  at  300  between  the  minimum  (50%)  and  the  maximum  (100%) 
values  allowable  for  the  variable. 
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Once  the  map  with  all  the  feasible  points  and  the  related  cost  is 
completed,  the  controller  determines  the  power  split,  in  particular 
the  power  to  be  sent  to  the  battery,  by  minimizing  the  cost  function 
Jdp  The  sample  time  for  the  computation  of  the  optimal  control 
policy  has  been  set  at  1  s  for  the  results  presented  in  this  article,  in 
order  to  achieve  a  low  computational  cost. 


4.  Simulation  results 

Figs.  17-19  show  the  simulation  results  during  the  Supple¬ 
mental  Federal  Test  Procedure  US06  (US06)  for  the  different  control 
algorithms  presented  in  Section  3.  The  RB  and  MPC  results  derive 
from  the  implementation  of  the  controllers  in  the  vehicle  simulator 


Time  [s] 


Fig.  19.  Battery  state  of  charge  during  the  US06  driving  cycle:  comparison  between  battery  only,  RB,  MPC  and  DP. 
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RB  MPC  DP 

■  NEDC  ■FTP  ■  ARTEMIS  ■  US06 


Fig.  21.  Percentage  of  life  increase  for  the  different  algorithms  and  driving  cycles. 


presented  in  Section  2.  In  order  to  fairly  compare  the  DP  algorithm 
results  with  those  from  the  other  controllers,  a  look-up  table 
describing  the  optimal  battery  power  profile  obtained  from  the  off¬ 
line  optimisation  procedure  has  been  derived  and  input  into  the 
simulation  model  of  Section  2.  The  supercapacitor  initial  state  of 
charge  is  set  at  80%. 

The  battery  current  profile  together  with  the  peak  and  the  RMS 
values  for  the  different  controllers  along  the  US06  driving  schedule 
are  reported  in  Fig.  17.  It  is  worth  noticing  the  difference  between 
the  DP  and  the  other  controllers  for  positive  battery  currents  (i.e. 
recharge  condition),  in  which  the  DP  recharges  the  battery,  whilst 
the  RB  and  the  MPC  assign  all  the  recharging  power  to  the  super¬ 
capacitor  stack. 

The  supercapacitor  state  of  charge  for  the  different  control  al¬ 
gorithms  is  reported  in  Fig.  18.  The  plot  shows  that  the  DP  algorithm 


and  the  RB  controller  maintain  a  higher  supercapacitor  state  of 
charge  for  a  large  portion  of  the  driving  schedule  compared  to  the 
MPC.  These  comments  can  be  inverted  for  the  battery  state  of 
charge  reported  in  Fig.  19,  for  which,  however,  the  DP  produces  a 
higher  final  state  of  charge  as  a  consequence  of  the  recharging 
power  sent  to  the  battery  stack,  as  discussed  earlier. 

Fig.  20  shows  the  normalised  battery  wear  for  the  different 
controllers  (RB,  MPC  and  DP)  along  the  same  driving  cycle.  As  ex¬ 
pected,  DP  represents  the  optimal  solution  achieving  the  lowest 
value  of  battery  ageing  (corresponding  to  Aheff,Dp)  despite  the  sig¬ 
nificant  battery  wear  profile  in  the  first  portion  of  the  driving  cycle. 
The  MPC  generates  a  lower  normalised  battery  wear  until  the 
supercapacitor  state  of  charge  approaches  its  lower  boundary. 
Starting  from  that  point,  the  MPC  battery  wear  significantly 
increases. 

The  predicted  battery  life  increase  is  summarised  in  Fig.  21. 
Table  6  reports  the  detailed  results  for  the  New  European  Driving 
Cycle  (NEDC)  and  for  the  Federal  Test  Procedure  (FTP75)  driving 
cycle,  while  Table  7  reports  the  overall  results  for  the  ‘Assessment 
and  Reliability  of  Transport  Emission  Models  and  Inventory  Sys¬ 
tems’  (ARTEMIS)  and  US06  driving  cycles.  The  previous  tables  are 
derived  for  nominal  thermal  conditions  of  the  drivetrains  and  the 
HESS.  The  adoption  of  the  supercapacitor  significantly  extends 
battery  life,  reducing  the  RMS  and  maximum  battery  current,  but  it 
does  not  improve  the  overall  efficiency  of  the  system  (notice  the 
values  of  power  loss  in  the  tables),  mainly  because  of  the  power 
losses  in  the  DC/DC  converter. 

The  supercapacitor-to-HESS  power  split  distribution  along  the 
US06  driving  cycle  is  shown  in  Fig.  22.  It  is  worth  noticing:  i)  the 
linear  distribution  between  the  supercapacitor  and  the  overall 
HESS  power,  suggesting  that  a  set  of  rules  might  be  derived  from 


Results  for  the  different  options  along  the  NEDC  and  FTP75  driving 
Parameter  NEDC 


Battery  only  RB 


Max  battery  current  120  106 

RMS  battery  current  41  14 

Battery  dissipation  205  24 

DC/DC  converter  dissipation  -  187 

Supercapacitor  dissipation  -  18 

Battery  life  increase  -  637 


schedules. 


MPC 


71 

25 

219 

20 

618 


DP 


23 

12 

17 

240 

25 

877 


FTP75 

Battery  only  RB 

127  106 

42  18 

331  60 


16 

378 


MPC  DP 


41  19 

10  10 

19  17 

445  477 

34  42 

1334  1481 


[Wh] 

[Wh] 

[Wh] 
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Fig.  22.  Power  split  distribution  along  the  US06  driving  cycle. 


the  DP  algorithm,  and  ii)  the  presence  of  points  with  a  positive 
supercapacitor  power  in  excess  of  the  corresponding  HESS  power, 
demonstrating  that  the  battery  supplies  extra  amounts  of  power  in 
order  to  recharge  the  supercapacitor,  as  confirmed  by  Fig.  23.  This 
condition  is  neither  discussed  nor  implemented  in  RB  strategies  or 
other  controllers,  which  send  the  supercapacitor  stack  only  the 
power  available  from  regenerative  phases. 

Linear  power  distributions  also  result  from  the  DP  algorithm 
applied  to  the  other  (i.e.  ARTEMIS,  FTP75  and  NEDC)  driving  cycles, 
here  omitted  for  brevity.  A  linear  approximation  of  the  scatter  plots 
of  the  power  split  distributions  can  be  expressed  as  in  Eq.  (46): 


pscap  =  aPHESS  +  b  (46) 

The  values  of  a  and  b  for  the  driving  cycles  simulated  in  this 
paper  are  reported  in  Table  8. 

5.  Conclusions 

Three  control  algorithms  for  the  power  split  within  a  HESS  of  a 
TTRP  HEV  have  been  presented,  two  of  them  being  newly  based  on 
a  battery  wear  model.  The  simulation  results  have  been  discussed 
and  analysed  to  assess  the  potential  benefits  deriving  from  the 


Fig.  23.  DP  power  split:  detail  of  the  US06  driving  schedule. 
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Table  8 

Linear  regression  of  the  power  split  between  HESS  and  supercapacitor. 


a H_ b  [W] 

NEDC  0.94  2376 

FTP  0.97  1997 

ARTEMIS  0.96  3657 

US06  0.97  2964 

Average  0.96  2749 


adoption  of  a  HESS.  In  particular,  the  study  demonstrates  that  a 
significant  decrease  of  battery  wear  and  RMS  and  peak  values  of 
battery  current  can  be  achieved  through  hybrid  energy  storage, 
whilst  the  adoption  of  the  supercapacitor  does  not  improve  the 
overall  energy  efficiency  of  the  system  in  nominal  thermal  condi¬ 
tions,  because  of  the  losses  in  the  DC/DC  converter.  The  rule-based 
controller  allows  a  67%  reduction  of  the  RMS  values  of  battery 
current  along  a  selection  of  driving  cycles  in  comparison  with  the 
same  vehicle  equipped  with  battery  only.  In  the  same  conditions 
the  battery  peak  current  is  reduced  by  38%.  The  model  predictive 
controller  and  the  dynamic  programming  algorithm  bring  an 
additional  reduction  of  the  root  mean  square  value  of  6%  and  10% 
respectively,  whilst  the  peak  values  are  additionally  decreased  by 
17%  and  45%. 

The  recommended  future  developments  are:  i)  to  define  a  set  of 
rules  starting  from  the  DP  algorithm  results  in  order  to  design  an 
enhanced  rule-based  controller,  and  ii)  to  implement  a  computa¬ 
tionally  efficient  DP-derived  algorithm  using  the  information  from 
the  cloud  and  the  vehicle  infotainment  system. 
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